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Abstract 



The non-Abelian discrete symmetry D-j of the heptagon is successfully applied to 
both quark and lepton mass matrices, including CP violation. 



1 Introduction 



The structure of quark and lepton mass matrices has been under theoretical study for many 
years. Whereas the 6 quark masses and the 3 mixing angles and 1 CP violating phase in 
the quark sector are now measured with some precision, the lepton sector is still missing 
some crucial information. Recently, the neutrino mixing angle #13 has been measured by 
the Daya Bay [1] and RENO [2] collaborations. The fact that sin 2 2# 13 is now centered at 
around 0.1 means that the previously favored tribimaximal mixing pattern (sin 2 #23 = 1/2, 
sin 2 ^ = 1/3, #i3 = 0) is invalid, although the A* symmetry [3JIUI2] used to obtain it [6] is 
still applicable with some simple modifications [3 13 [9]. On the other hand, in the simplest 
application [3j [5] of A4, all the quark mixing angles are zero. The question is whether there 
exists another symmetry which successfully yields both quark and lepton mass matrices, 
with good fits of all masses, mixing angles, and phases. The answer is yes, as elaborated 
below. 

Using the non-Abelian discrete symmetry Dj of the heptagon, it has been shown p.0] that 
the CP violating phase of the quark mixing matrix may be predicted, whereas D-j also yields 
a pattern [11] for the neutrino mass matrix consistent with what is observed. This pattern 
is previously derived using the symmetry Q$ [12], and realizes a specific conjecture [13] that 
the neutrino mass matrix has two texture zeros in the basis that charged-lepton masses are 
diagonal. 

In Sec. 2 the symmetry D 7 is explained. In Sec. 3 the assignments of quarks under 
D 7 are given with the accompanying Higgs structure and the resulting mass matrices. In 
Sec. 4 numerical fits to the quark masses and mixing angles are given, with a prediction of 
the CP violating phase. In Sec. 5 the assignments of leptons under D 7 are given with the 
accompaning Higgs structure and the resulting mass matrices. In Sec. 6 the neutrino mass 
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matrix is analyzed to show that it allows for nonzero #13 and a specific correlation between 
it and #23 as well as Sep- Given that 9 12 is close to the tribimaximal value, it prefers an 
inverted hierarchy of neutrino masses although a quasidegenerate pattern with either normal 
or inverted ordering cannot be ruled out. In Sec. 7 there are some concluding remarks. 

2 Heptagonic Symmetry D7 

The group D-j is the symmetry group of the regular heptagon with 14 elements, 5 equivalence 
classes, and 5 irreducible representations. Its character table is shown below. 



class 
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Xi 


X2 


X3 


Xi 


X5 


Ci 


1 


1 


1 


1 


2 


2 


2 


c 2 


7 


2 


-1 


1 











c 3 


2 


7 


1 


1 


d\ 


a 2 


a 3 


c A 


2 


7 


1 


1 


a 2 


a 3 


ai 


c 5 


2 


7 


1 


1 


a 3 




a 2 



Table 1: Character Table of D-j. 



Here n is the number of elements and h is the order of each element. The numbers 
are given by au = 2 cos(2A;7r/7). The character of each representation is its trace and must 
satisfy the following two orthogonality conditions: 

n iXaiX*bi = n5 ab , n iXaiX*aj = TlSij, (1) 

Ci Xa 

where n = J2i n i is the total number of elements. The number of irreducible representations 
must be equal to the number of eqivalence classes. 

The three irreducible two-dimensional reprsentations of D 7 may be chosen as follows. For 
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2i, let 

Ci : (J J), C 2 : (J_ fc (/c = 0, 1,2,3,4,5,6), 

/ w 0\ / w « 0\ /u, 2 0\ 0\ 

-to °mi 

where w = exp(27ii/7), then 2 2) 3 are simply obtained by the cyclic permutation of 6*3,4,5. 

For D n with n prime, there are 2n elements divided into (n + 3)/2 eqivalence classes: C\ 
contains just the identity, C 2 has the n reflections, from k = 3 to (n + 3)/2 has 2 elements 
each of order n. There are 2 one-dimensional representations and (n — l)/2 two-dimensional 
ones. 

The group multiplication rules of L> 7 are: 

1' x 1' = 1, 1' x 2, = 2,, (3) 
2, x 2, = 1 + T + 2 m , 2, x 2 m = 2, + 2 i+2 , (4) 

where 2 4j5 means 2i j2 . In particular, let (ai,a 2 ), (61,62) ~ 2i, then 

ai6 2 + a 2 6i~l, ai& 2 - a 2 h ~ 1', (ai&i, a 2 6 2 ) ~ 2 2 . (5) 

In the decomposition of 2i x 2 2 , we have instead 

(a 2 b 1 ,a 1 b 2 ) ~ 2i, (a 2 &2,ai&i) ~ 2 3 . (6) 

3 Quark Sector 

We assign quarks as shown in Table 2 and Higgs doublets as shown in Table 3, together with 
an extra Z 2 x Z 2 symmetry. 
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symmetry 


(u,d), (c,s)\ 


(t,b) 


(d c , s c ) 




(w c , c c ) 


t c 


D 7 


2i 


1 


2i 


1 


2 2 


1 


zi 


+ 


+ 






+ 


+ 




+ 


+ 


+ 


+ 




+ 



Table 2: Quark assignments under D 7 x Z$ x Z%. 



symmetry 
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$ 2 


$3,4 


$5,6 


$7,8 


D 7 


1 


1 


2i 


2 2 


2 3 


zi 


+ 






+ 


+ 


2% 


+ 


+ 


+ 


+ 





Table 3: Higgs doublet assignments under D 7 x Z$ x Z%. 



As a result, the (u, c, t) mass matrix is diagonal, coming from the Yukawa terms uu c <j) 7 + 
cc c 0g and tt c ^{. As for the (d, s, b) mass matrix, the allowed Yukawa terms are (ds c + sd c )4>2, 
bb c <f)2, b(d c (f)\ + s c (f)^), and {dcf)\ + s<f)^)b c . The resulting mass matrix is thus of the form pU] 

/ a f&\ 

a b , 
\fc c d / 



M d 



(7) 



where 



4 Prediction of CP Phase 



As in Ref. [10], we can redefine the phases of M.d so that a, b, c, d are real, but £ is complex. 
Since Ai u is diagonal, we have 



/m d \ 

m s 
V mj 



\ 



m? 



\ ml) 
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where Vl is the observed quark mixing matrix up to phase conventions. The structure of 
Aid-Mj allows us to obtain the following first approximations: 

ac + £,bd 



where a 2 << b 2 and |£| 2 << 1 are assumed. We now rotate Aid-M^, using 



V* 



I 1 o v ub \ 

1 v cb 
\-v: b -v; b i ) 



to obtain the 2x2 matrix 



where 



M 2 M\ 



A C\ 
C* B ) ' 



A 
B 
C 



yielding 



IK. 



= a 2 + \i\ 2 b 2 -\V ub \ 2 m 2 b , 
= a 2 + b 2 -\V cb \ 2 ml 
= ib 2 -V ub V: h ml 

l(B + A) + y(B-Ay + 4\C\ 2 , 
1 1 



1 - 



4|CP 



2 2^ (B - A) 2 + A\C\ 2 ' 
where the phase of V us is that of C, and 

m d = \2abc£ — a 2 d\/m s m b . 



(9) 



(10) 



(12) 
(13) 
(14) 



(15) 
(16) 



(17) 



Using \V US \ = 0.22534, we find \C\ 2 /(B - A) 2 = 0.05971, and m 2 » m 2 implies A ~ 
0.053515, hence m 2 ~ 1.053495. Using these formulas, the 6 parameters a, 6, c, d, Re(£), im(£) 
may then be determined and the CP violating parameter J is predicted. 
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Figure 1: The CP violating parameter J versus m^/m^. The solid (dash) lines indicate the 
one (two) standard-deviation bounds of J. 

For our numerical analysis, we dispense with the approximate expressions and diagonalize 
Aid-M-li directly. We scan for solutions consistent with data on the 3 masses and 3 mixing 
angles, within one standard deviation of each parameter. We then obtain J numerically from 
the resulting Vckm- This is then the prediction of our model. In Fig. 1 we plot J versus 
m s /md, which shows good agreement with data. We use the 2008 updated values [H] of 
m d,s,b evaluated at M w : 

m d (M w ) = 2.93 (+1.25/- 1.21) MeV, (18) 
m s (M w ) = 56 ±16 MeV, (19) 
m b (M w ) = 2.92 ± 0.09 GeV, (20) 

and the 2012 Particle Data Group (PDG) [15] values of the mixing angles 

\V US \ = 0.22534 ±0.00065, (21) 
\V cb \ = 0.0412 (±0.0011/ - 0.0005), (22) 
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\V ub \ = 0.00351 (+0.00015/ -0.00014). (23) 

Note that PDG also lists the condition 17 < m s /m d < 22 and the value of the CP violating 
parameter is 

J = 2.96 (+0.20/ - 0.16) x 10" 5 . (24) 

We show in Table 4 sample values of a, b, c, d, Re(£), Im(£) with the corresponding values of 
m d ,m s ,m b , \V US \, \V ub \, \V cb \ and J as well as m s /m d . 



a (GeV) 


b (GeV) 


c (GeV) 


d (GeV) 




MO 


m s /m d 


m d (MeV) 


m s (MeV) 


m b (GeV) 


\v us \ 


\v ub \ 


\Vcb\ 


J 


0.0125 


0.138 


1.32 


-2.60 


0.053 


-0.084 


17.00 


3.89 


66.2 


2.92 


0.22534 


0.00355 


0.0420 


2.95 x 10~ 5 


0.0124 


0.139 


1.34 


-2.60 


0.058 


-0.084 


17.25 


3.91 


67.4 


2.93 


0.22532 


0.00358 


0.0420 


2.89 x 10~ 5 


0.0123 


0.138 


1.40 


-2.60 


0.064 


-0.087 


17.50 


3.96 


69.2 


2.96 


0.22519 


0.00363 


0.0409 


2.76 x 10~ 5 


0.0122 


0.138 


1.39 


-2.55 


0.068 


-0.084 


17.75 


3.94 


69.9 


2.91 


0.22501 


0.00359 


0.0415 


2.70 x 10~ 5 



Table 4: D 7 parameter fits of quark masses and mixing. 



5 Lepton Sector 

Using again D 7 x Z$x Z%, we assign leptons as shown in Table 4 and Higgs triplets as shown 
in Table 5. 
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symmetry 


(v e , e) 




e c 


[(/Ar c )] 


D 7 


1 


2i 


1 


2 3 


Z* 


+ 


+ 


+ 


+ 


z% 


+ 


+ 


+ 


+ 



Table 5: Lepton assignments under D 7 x Z 2 d x Z 2 " . 



symmetry 


6 


£2,3 


D 7 


l 


2i 


Z* 


+ 


+ 


z$ 


+ 


+ 



Table 6: Higgs triplet assignments under D 7 x Z% x Zg. 



As a result, the (e, fi, r) mass matrix is diagonal, coming from the Yukawa terms ee c 0° 
and /x/x c 05 + TT c 0g. As for the Majorana (i/ e , u^, u T ) mass matrix, the allowed Yukawa terms 
are v e vj^[, (z^z/ T + v T v^)^[, and ^(^£3 + ^V*^)- The resulting mass matrix is thus of the 
form 111 



(a c d\ 
M u = c b , 
K d b 0) 

which was first derived using Q 8 [12] . and realizes one of the conjectures of Ref. [13]. 



(25) 



6 Analysis of Neutrino Mass Matrix 



Rotating M. v to the tribimaximal basis using 



it becomes 
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— u TB 
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\A/3 






W) 
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( m i 








M {1,2,3) _ 




m 6 


m 2 m 5 














m 5 m 3 / 
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(26) 



(27) 



where 



rrii 
m 2 
m 3 

m 5 
m 6 



(2a + b- 2c- 2d), 
(a + 2b + 2c + 2d), 



-c + d), 



3 

-6, 
1 

Te [ - c + d) = T2' 

1 -.(2a - 2b + c + d) = —=(mi + 2m 2 + 3m 3 ). 



(28) 
(29) 
(30) 
(31) 

(32) 

(33) 



3^2 2^2 
If 777.4 = = W6 = 0, tribimaximal mixing is recovered. In particular, 7774 7^ or 7715 7^ 
means that 6 Vi 7^ 0. In previous studies, the special cases m 4 7^ 0,m 5 = m 6 = [HI [17] 
and 7775 7^ 0,7774 = 7776 — [H [91 [18] have been explored. The requirement from D7 that 
7775 = 7ti 4 /a/2 is a new condition which will predict a special correlation between 13 and #23 
as well as <5cp. 



Consider the unitary matrix U e such that 

^(1,2,3)^(1,2,3))^ = 



/ 


777^ 2 








\ 







1 / 12 

|m 2 | 2 







V 








1 777.3 


V 



(34) 



then U' ai = UtbU 6 is the lepton mixing matrix up to phases. Let U e be approximately given 
by 



U r 



then for |m^| ~ |mi| , we have 



£21 



/ 1 e i2 ei 3 \ 

^21 1 e 23 

\e3i £32 1 / 



— (7776777^ + 7772777g) 
|7772| 2 — |777i| 2 



(35) 



(36) 



In addition, since the effective neutrino mass m ee in neutrinoless double beta decay is given 
by 



777 ee = \a\ = |777i + 777 2 + 777 3 | 
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(37) 



whereas 



we have the relationship 



m 3 = ]-(2V2m 6 - mi - 2m 2 ) , (38) 



\m 3 \ 2 -m 2 ee = ^(|m 2 | 2 - |mi| 2 )[l + 4y/2Re(e 21 )]. (39) 

Since |m 2 | 2 — \mi\ 2 ~ Am^ is very small, this model predicts m ee = \m 3 \ to a very good 
approximation. The structure of Eq. (38) also shows that an inverted ordering of neutrino 
masses is expected, although the quasidegenerate limit is also possible in which case either 
inverted or normal ordering may occur. In the following we focus on the inverted case, i.e. 

1^3 1 < \ m l\ < |?™2| ■ 

For m 4 7^ 0, v 3 is rotated to u' 3 according to 

m\m\ + m 4 mn m 2 m! + m^m\ . , 

e i3 - 19 lz -9. » e 23 - /k„„, ,, - , 2V ( 40 ) 



|m 3 | 2 -|mi| 2 ' "" y/2(\m 3 \ 2 - \rrix 



As a result, 



/2 /T -m 4 (mi + 2m 2 )* +m*(2 mi + m 2 ) 

- \/3 £ - + \/3 £ -" V5(M»-KI») ' (41) 



l 1 1 1 (mi — m 2 )m 4 

°" = "71 " + vf 23 = "2 ~ VB(| mi p - kF)' (42) 



r ,, 11 1 1 [m x -m 2 )m\ 

u « " 7S"7g 6l3 + 75 e ' 23 "71" V6(| mi p-Kl') ' (43) 

If all parameters are real, then for U' e3 = 0.16, sin 2 29 23 would be 0.80, which is ruled out by 
present data, i.e. sin 2 2 9 23 > 0.92. However, a fit may be obtained for complex values. 

We go back to Eq. (25) and observe that a, c, d may be chosen real, so only b is complex. 
This means that m 4 is real as well as 2m 1 — m 2 , and for m 6 = 0, m 3 = — (mi + 2m 2 )/3. 
Writing mi j2 as mi j2 e^ 1 - 2 with m 2 ~ mi and sin0 2 = 2sin0i, we obtain 

U' e3 ~ —1=- — — [- cos 0i + cos0 2 - 9«sin0i], (44) 



V^Am 2 



32 

11 



1 



17111714 



im| 2 



[cos 0i — cos 02 — i sin 0i] 



rr' 1 mim4 r jl . • • ; -i 

U t3 ~ - 7 = + —, =t — - cos 0i -cos02 -zsm0ij, 

V2 v6Am3 2 



(45) 
(46) 



where cos 02 = ± ^/l — 4 sin 2 0i. We then have 



sin 2 #13 = 





U' e3 


2 


1 + 


ei3 


2 + 


^23 


2 



TV I 2 



tan 2 #23 = ^. 



Since 



\mj\ ~ 



Am| 2V /5 + 4cos(0 2 - 0i) 



2yl - cos(0 2 - 0i) 

the above equations relate |m3| = m ee with #i 3 and # 2 3- If we fix #13, we then obtain 
function of #23- We plot in Fig. 2 our model predictions for |mi j2 | and |ms| = m ee 
The other data points are taken to be their experimental central values. 



a 

sin 2 2# 23 



(47) 

(48) 

m 3 | as 
versus 



0.10 
0.08 



% 0.06 



0.04 



0.02 




0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 

sin 2 (2 # 23 ) 

Figure 2: Neutrino masses m\ 2 and m 3 = m ee versus sin 2 26* 2 3- 
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If we rotate M]; 2; \M]; 2 ^ by 



IK 



( 1 e 13 \ 





V- £ 13 



^23 



(49) 



^23 1 / 

we obtain the 2x2 mass-squared matrix spanning u[ 2 . This differs from the 2x2 submatrix 
in the tribimaximal basis by terms quadratic in which are important in obtaining the 
correct Am^ and Eq. (36) becomes modified. However, we can adjust \m2\ versus |mi| as 
well as rriQ to fit the data. These adjustments will have negligible effects on \m 3 \. 

1.00 




0.93 0.94 0.95 0.96 0.97 0.98 0.99 

sin 2 (2 6 23 ) 

Figure 3: The CP violating parameter | sin<5cp| versus sin 2 26*23- 

We plot in Fig. 3 our model prediction for | sin5cp\ versus sin 2 20 2 3- To obtain smScp, 
we use 

U' 



1 



from which we find (using U' 2 = I ^"^2 1 e ie * M2 , etc. 



cos 6*23 sin 6 ~ \U' 2 \ sm(0 M3 - 6» m2 - e3 ). 



(50) 



(51) 
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7 Concluding Remarks 



We have studied a specific pattern for both quark and lepton mass matrices. In both cases, 
one mass matrix is diagonal (Ai u and Ai e ), whereas the other has two zeros (A4d and M. v ). 
In the case of Ai u , the assumption that it is Majorana corresponds to one of the conjectures 
of Ref. [13j, whereas the Dirac mass matrix hA d requires further restrictions to make it 
predictive, as first proposed in Ref. [10] using the non-Abelian discrete symmetry D 7 . The 
conjectured form of M. v was first derived [12] using Q 8 , but it may also be obtained [H] 
using D 5 or D 7 . Here we consider as the unifying symmetry for both quarks and leptons. 

The CP violating parameter J in the quark sector is constrained in this model by 
md,m s ,m b , \V US \, \V u b\, \V c b\. Within one standard deviation of all six measurements, we ob- 
tain J in agreement with data. In the neutrino sector, we obtain |mi 5 2| as well as \m^\ = m ee 
as functions of sin 2 2#23 and also predict sin Sep as a function of sin 2 26*23- 
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